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Abstract. We investigate the diffusion of a grain boundary in a crystalline material. 
We consider in particular the case of a regularly spaced low-angle grain boundary 
schematized as an array of dislocations that interact with each other through long- 
range stress fields and with the crystalline Peierls-Nabarro potential. The methodology 
employed to analyze the dynamics of the center of mass of the grain boundary and 
its spatio-temporal fluctuations is based on over-damped Langevin equations. The 
generality and the efficiency of this technique is proved by the agreement with molecular 
dynamics simulations. 
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1. Introduction 

Understanding interface kinetics in materials is an important theoretical and practical 
problem, since this process influences the microstructure, such as the grain size, the 
texture, and the interface type. From the theoretical point of view, the study of 
processes that involve surface and interface properties gained significant interest in 
non-equilibrium statistical mechanics PQ[2]. In particular, dislocations [3, 4] and grain 
boundaries [5j E] provide a concrete example of driven elastic manifolds in random 
media [TJ. Other example of this general problem are domain walls in ferromagnets 
0, [9], flux line in type II superconductors pm fTTj . contact lines p2l 03] and crack 
fronts p^l [T5] . From the point of view of applications, understanding grain boundary 
kinetics has a great importance for polycrystalline materials, since the resulting grain 
microstructure determines material properties such as strength, hardness, resistance to 
corrosion, conductivity etc. [IB] . Hence the ambitious goal of these studies is to be able 
to control the microstructural properties of polycrystals. 

Several approaches have been employed in the literature to study grain boundary 
kinetics. Ref. [TTJ employs molecular dynamics (MD) simulations with appropriate 
interatomic interactions to study the diffusion of grain boundaries at the atomic scale 
[17] . The method allows to quantify the mobility of grain boundaries and to compare the 
results with experiments [IT] . While MD simulations provide a very accurate description 
of the dynamics, the method suffers from numerical limitations and it is difficult to reach 
the asymptotic regime. An alternative method is provided by the Langevin approach in 
which the grain boundary is assumed to evolve stochastically in an external potential 
[18]. The dynamics of the underlying crystalline medium enters in the problem only 
through the noise term (due to lattice vibrations) and the periodic potential (Peierls- 
Nabarro). Hence, the equations of motion of the atoms or molecules are not directly 
relevant. Indeed there is experimental evidence in supporting of separation of time scales 
in plastic flow [19J and it is thus possible to integrate out the fast degrees of freedom 
(atomic vibrations) and consider only the slow ones (dislocations position). 

Here we study the evolution of a grain boundary (GB) in a crystalline material 
by the Langevin approach. The GB is treated as an array of interacting dislocations 
performing a thermally activated motion in a periodic (Peierls-Nabarro) potential. 
Similar models have been employed in the past to study the conductivity of superionic 
conductors [20l[2T], the relaxational dynamics of rotators [22] and Josephson tunneling 
junctions [IB] . Notice that the crucial role played by long-range stresses is often 
disregarded in analyzing GB deformation. On the other hand, it has been shown in 
Ref. [5] that a surface tension approximation for the GB stiffness is inappropriate and 
one has to consider explicitly non-local interactions. The present model incorporates 
this effect in the equations of motion. 

We simulate the set of Langevin equations numerically to describe the GB 
kinetics and its fluctuations. The results are in good agreement with MD simulations 
[IT] and allow to clarify the origin of the short time deviations from the diffusive 
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behavior observed in Ref. [17J. In addition, a linearized version of the model can be 
treated analytically and the asymptotic results are found in good agreement with the 
simulations. The manuscript is organized as follows: in Sec. II we introduce the model, 
which is first studied in the flat GB limit in Sec. III. Sec. IV presents numerical 
simulations of the full flexible GB problems and Sec. V discusses the continuum theory. 
Sec. VI is devoted to conclusions. 



2. The model 



To study the GB dynamics we consider a phenomenological mesoscopic approach. We 
consider in particular the case of a regularly spaced low-angle grain boundary schema- 
tized as an array of straight dislocations that interact with each other through long- 
range stress fields and with the crystalline Peierls-Nabarro (PN) potential. The GB is 
composed by N dislocations where configurations are repeated ad infinitum because of 
periodic boundary conditions along the y direction. Each dislocation has Burger vector 
of modulus b parallel to the x axis and the distance between two adjacent dislocations 
along the y direction is fixed to be a. Each straight dislocation interacts with the lattice 
and with others dislocations through long-range stress fields. The effect of the lattice 
over each n-th dislocation can be decomposed as the sum of three contributions: 

• FpN(x n ) = ~- / ^ 2nr ^( ^fr*™ )' the force where A is the area of the GB, /x is 
the shear modulus and r$ the inter-atomic distance; 

• -7x n (t), the average effect of the lattice fluctuations where 7 is the viscosity 
coefficient; 

• jr] n (t), the impulsive effect of the lattice fluctuations assumed to be Gaussian 
for the central limit theorem and uncorrelated in space and time: (i]n{t)) = 0, 
(r) n (t)r) m (f)) = DS nm S(t — t') where D is the diffusion coefficient [IS] . 



The long-range stress field exercised by all the other dislocations over the n-th, the 
Peach-Koehler force Fp^(x, y), is computed considering the image dislocations method 
to comply with periodic boundary conditions along the y direction. Making use of 
calculations in [231 121] one can find the following expression 



F n > N (~ y\ - „ rf> E ^ N (r -r V 

(1) 

{cosh[27r(x w - x m )/Na] cos[2Tc(y n - y m )/Na] - 1} 
{cosh[27r(x n - x m )/Na] - cos[2?r(?/„ - y m )/Na]} 2 ' 

where u is the Poisson's ratio, y n = n ■ a and y m = m ■ a. Finally the over-damped 
Langevin equation [IS] for the GB reads 

7i n (t) = F PN (x n ) + Fp'x (x, y) + yr} n (t), (2) 
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for n = l N, or rather 



2, 



x n (t) - -A^^ sin(^^) - jV 2 a 2 i ^ 1 1 ^ 7 E™=i(^ - 

{cosh[27r(x n - x m )/Na} cos[27r(?/ n - y m )/Na] - 1} , v 
' {cosh[27r(x n - x m )/iVa] - cos[27r(y n - y m )/Na}} 2 ^ ln{ h 



(3) 



To indicate the amplitude of the Fpn and Fpk forces we introduce respectively the 
parameters A PN = Aiibjlur^ and A PK = fib 2 7i/a 2 (l — u)j. 

The key quantities that we consider in order to characterize the dynamics of the 
GB are: 

• the mean-square displacement of the center of mass, Ax cm (t) = (x 2 m (t)) — 
(x cm (t)) 2 = (x n (t) ) - (x n (t)) 2 where x cm (t) =x^= l/iV^ =1 i„(t); 

• the mean-square width W 2 {t) = (x^) — ((x^) 2 ). 

In the following, we first analyze the case of a flat GB for which a comparison with 
MD simulations approach [17] is made. Next we consider the full flexible description 
of the GB. Finally, we discuss a linearized version of the model that can be treated 
analytically. 



3. Flat grain boundary 

For many applications a good approximation is to consider a flat GB with a single 
degree of freedom, for which Fp^(x, y) = and x n (t) = x cm (t) for n = 1, N. In 
other words, the flat GB is described by the following equation 

x cm (t) = sin(^L) + 77(i), (4) 

where the correlation properties of the thermal fluctuations are: (^(t)) = and 
{r](t)r](t')) = D8(t—t'). This type of equation, also known as the Kramers equation with 
periodic potential, has been extensively studied in the literature [18]. In particular, the 
mean-square displacement is known to display a combination of oscillatory and diffusive 
behavior [TS1122]. Different dynamical regimes are found as the potential strength or the 
friction varies [25]. In fact, we show next that this simple model allows to understand 
the short-time deviations from diffusive behavior observed in MD simulations [T7] . 

Integrating Eq. |4] with the initial condition x cm (0) = by means of computer 
simulations (fitting D and 7 with the condition Dj = Dr/A4 [H] where M. is 
the mobility and Dr is the renormalized diffusion coefficient considered in [17]), we 
have compared the mean-square displacement Ax cm (t) to the one obtained from MD 
simulation in Ref. [17] . In Fig. [1] this comparison is displayed together with the 
mean-square displacement of the renormalized free Brownian motion (described by the 
equation: x = rj(t) with (i](t)) = and (i](t)r](t')) = Dn5(t — t')). The agreement 
between the two simulations is extremely good. For higher times (t > 80ps), the 
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mean-square displacement tends to the renormalized Brownian motion. Hence taking 
explicitly into account the sinusoidal Peierls-Nabarro force in the Langevin equation 
allows to describe the mean-square displacement for early times of the dynamics. 
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Figure 1. Mean-square displacement Ax cm (t) of the fiat grain boundary comparison 
between molecular dynamic simulation |17j (green line) and Langevin approach 
simulation (dashed blue line). Both type of simulations predicts a linear dependence 
in time of Ax cm (t) for long times represented in the figure by the renormalized free 
Brownian motion (straight black line). 



One can deduce in a simplified intuitive way the temporal evolution of the mean- 
square displacement starting by the transition probability density P for small times 
(small r) [18] 

Y [x-x'-F pn (x)t] 2 
P(x,t + T\x',t) = ; e 4Cr . (5) 

Next we consider the transition probability Poi(x + Axi,t + r\x ,t ) to run from 
the point Xq at time to to the point X\ = Xq + Ax\ at time t\ = to + r and 
Pi2(%o + Ax% + Ax2, to + 2t\xq + Ax±, to + r) to run from x\ at t\ to X2 = x\ + Ax2 at 



U +T 



P 



01 



p 



1 
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[Axi-F pn (x )t} 
ADt 



12 



[Ax 2 -F pn {x q +Ax 1 )t} 2 



(6) 



2 V tcDt 

For a free Brownian motion (Fpj^(x) = 0) the condition P i = Pu implies Ax% = Ax2 
and stochastic displacements are space independent. If we impose this condition in 
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presence of a periodic force F PN (x) = —dUp^ix) / 'dx, we obtain 

Poi = P12 Ax 1 - F pn (x q )t = Ax 2 - F PN (x + Axi)r 



and then 



Ax2 = Ax\ 



Ax2 ^ Ax\ if 



1 + 



dF PN (x) 



dx 



T 



.T() 



(7) 



dF, 



PN 



dx 



^ or rather 



d 2 U 



PN 



dx 2 



£0. 



This result implies that, with the initial condition x cm (0) = 0, if the potential Up^(x) is 
convex (concave) the mean-square displacement curve is concave (convex). In the case of 
the PN potential, we find indeed that the mean-square displacement curve should display 
upper and lower deviations from the straight line, corresponding to a renormalized free 
Brownian motion, depending in d 2 UpN(x)/dx 2 . These deviations decrease with time so 
that for large times the curve should approach a straight line |18j . 



4. Flexible grain boundary 

A more general description of the GB considers its internal deformation and the 
dynamics is described by Eq. [31 The dynamical behavior of the GB depends on the 
amplitude of the three terms in the right-hand side of Eq. [31 The parameters that 
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Figure 2. Mean-square width of the grain boundary, W 2 (t), in the case of N = 32 for 
the two typical situation in which (continuum line) the GB exfoliate (W 2 (t) increase 
with time) because the noise is high enough respect to Apk and to Apn, (dashed line) 
the GB reach a stationary state (W 2 (t) saturates after a certain time) because the 
noise is small enough respect to Apk or to Ap^- 



characterize the behavior of the GB are a, 6, Apk, Apn and D. Varying the values 
of these parameters, in the long-time limit the GB can either exfoliate (when the 
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noise, D, is high enough with respect to A PK and to A PN ) or reach a stationary state 
(when the noise, D, is small when compared to Apx or to Ap^). The asymptotic 
behavior can be read off from the width W 2 (t) that keeps on increasing when the 
GB exfoliates and saturates when the GB remains stable. In Fig. [2] the comparison 
between these two typical situations is displayed in the case of N = 32, a = b = 3n, 
A PN = 0, 0.1, 0.2, 0.4, D = 0.25 and A PK = 0.14 for the case in which the GB remains 
stable, while Apx = 0.0896 for the case in which the GB exfoliates. 

In what follows, we analyze the dynamical behavior of the stable GB for a = b = 3n, 
A PN = 0,0.4, A PK = 0.14 and D = 0.25. In Fig. [3] the average position of the GB 
center of mass Ax cm (t) is displayed with and without the PN force in Log-Log scale for 
N = 32,64,128,256,512. 




(b) 



Figure 3. Mean-square displacement of the center of mass of the grain boundary, 
Axcm(t) , for Langevin approach simulation without the PN force (a) and with the 
PN force (b). Ax cm {t) for N = 32, 64, 128, 256, 512 is displayed in Log-Log scale. 



For long times in both cases we have a linear behavior Ax cm (t) ~ t, but for short 
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Figure 4. Mean-square width of the grain boundary, W 2 (t), for Langevin approach 
simulation without the PN force (a) and with the PN force (b). W 2 (t) for N = 
32, 64, 128, 256, 512 is displayed in Log-Log scale. 

times, in the presence of the PN force, there is a clear deviation from linearity. This 
result confirms the conclusion made in the previous section, that the PN force is the 
cause for the deviation from linearity of Ax cm (t) for short times observed in Ref. |17j . 
Next we characterize the morphology of the GB through the width W 2 (t). In Fig. 0] 
W 2 (t) for N = 32, 64, 128, 256, 512 is displayed with and without the PN force in Log- 
Log scale. In the absence of the PN force (Fig. 0^.) the time dependence of W 2 (t) is 
qualitatively similar to the same case but with linearized PK force discussed in the next 
section, while in the presence of the PN force, for Ap^ = 0.4 (Fig. 0b), W 2 (t) exhibit a 
plateau for intermediate times. 
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5. Continuum Theory 

It is possible to develop an analytic expression in the continuum limit (a — > 0, N — > oo 
and L = iVa = const.) for short or long times for W 2 (t) in absence of the PN force 
linearizing the PK force. The equation of motion for Fpn = and Fpx linearized is 

in(t) = -7T7T- — ^E?"^ — -^+Vn(t). (9) 
2tt(1 - z/) 7 ^ (y n - y m y 



To obtain the short time behavior is sufficient to rewrite Eq. [9] as a generalized 
Ornstein-Uhlenbeck process (18] 

N 

Xn{t) = E 9nmXm + Vn(t) . (10) 
m=l 

The general solution of Eq. [10] is 

X n(t)= / E Gnm{t')T] m {t - t')dt' , (11) 

with {G nm {t)} = G(t) = e §t = 1+gt+gH 2 /2 + ... (where 1 = {<%}). From the definition 
of W 2 (t), results 

n JL rt D N r* 

w2 w = ^ E / - E / GUVGUW- (12) 

Replacing the Taylor expansion of the G matrix in Eq. [12] one obtains for short 
times (t < 1/II^H) that iy 2 (t) = (1 - l/N)Dt + o(t) and in the continuum limit 
W 2 (t) = Dt + o(t). To obtain the long times behavior of W 2 (t), we rewrite Eq. [9] 
in Fourier space [5] [26]. Employing the decomposition x m = 1/ L^2 k exp(—ikam)xk, we 
obtain 

,,1,2 °° °° ( —ik'am —ik'an\~, 

= - o n ^ 2f E e " am E Lfc ' ( 7 v, } * + (13) 

27r(l — z/WgtL ^— ' (m — n) z 

m=— oo n=— oo ' 

The first term in the right-hand side of Eq. [13] can be rewritten as 

Z 10 \ „ \ A(k-k')amt \ 1 e 



2tt(1 - z/) 7 a 2 L ^ ^ v ^ d 2 

v ' fc' m=— oo <i=— oo 

where d = m — n. Using the following results 



El 7T COS{Cd) 7T 7T|C| C 

we obtain 



d 2 6 ' ^ d 2 6 2 4 

d=l d=l 



~ l_ e ^ o ^ 1 - cos(fc'ad) „„ A;' 2 a 2 
E ^2 =2 E ^ = 7r|fc|a-— , (16) 



d? d 2 

d=—oo d=l 



so that Eq. [13] becomes 



ik = -2^r^ m ~ ~ )Xk + Vk ' (17) 
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In the long time limit (large x, small k) the k 2 term can be neglected. Finally 
in the continuum limit we replace y n , y m by the continuum variables y, y' and 
( V (y, t) V (y', f)) = aD5(y - y')5(t - t'). 
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Figure 5. Comparison between the mean-square width W 2 (t) for linearized Fpk and 
Fp 7v — in the case of N = 32 computed by numerical simulation and the continuum 
theoretical prediction in the short and long times limit. The black line represents 
the simulated data, the green line the long times theoretical prediction with fitted 
parameters and the blue line the short times theoretical prediction. 
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Figure 6. Size dependence for the saturation value of the mean-square width W 2 
in Log scale for the abscissa (x-axis). As can be seen in figure, in the case in which 
the PN force is absent {Apm = 0) the relation is logarithmic (W 2 ~ logiV), while in 
presence of the PN force (Ap^ 7^ 0) a slight deviation from the logarithmic dependence 
is observed. 
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Thus Eq. H3] becomes 

x k = -a\k\x k + rj k , (18) 
where a = /x6 2 /[2(l — u)^a 2 }. Eq. [TBI can be solved exactly and W 2 {f) is given by [21 127] 

w2 = aD_n,L + j n(1 _ g-^rot/LNl (19) 

27ra a 

To reproduce the continuum limit by simulations of Eq.[9]we would need a very large 
GB, with N 3> 512. Thus a comparison between Eq. [19] and the simulations for small 
N is possible only by introducing some effective parameters in Eq. [19] In Fig. [5] W 2 (t) 
computed by the simulations with N = 32 (for which the better statistic is available) 
is compared with the fitted theoretical prediction for short and long times. Eq. [19] also 
predicts that the saturated value of width (W 2 ) exhibits a logarithmic dependence on the 
GB length L = Nf a: W 2 ~ log N. In the case in which Fp^ is not linearized this result 
is also confirmed by numerical simulations, showing that W 2 increases logarithmically 
with N when F PN = (see Fig. [6]). In presence of a periodic potential (F PN > 0), 
however, we observe a deviation from the logarithmic growth at large N. This suggests 
that the Peierls-Nabarro potential may set a limit to the GB roughness. 

6. Summary and Discussion 

We have investigated the diffusion of a regularly spaced low-angle grain boundary in 
a crystalline material. A typical computational method to describe the dynamics of 
the grain boundary is to perform deterministic molecular dynamics simulations with 
appropriate interatomic interactions [17J. Here we have employed the over-damped 
Langevin approach to obtain a long time description of the dynamics, but in particular 
to perform a comparison with molecular dynamics simulations for a specific material 
[T7] . The first results is the interpretation of the early times behavior of the mean- 
square displacement Ax cm (t). The deviation for early times of Ax cm (t) by the case 
of the renormalized Brownian motion, that holds for long times, can be interpreted 
as the effect on the dislocations of the periodicity of the lattice giving rise to the 
Peierls-Nabarro potential. Secondly the description of the dynamic (Ax cm (t)) and the 
morphology (W 2 (t)) of the grain boundary by means of over-damped Langevin equations 
is in qualitatively good agreement with its behavior in real materials, so this approach 
can be considered an useful tool for these studies. 
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